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ABSTRACT 

There are now about fifty known radio pulsars in binary systems, including 
at least five in double neutron star binaries. In some cases, the stellar masses 
can be directly determined from measurements of relativistic orbital effects. In 
others, only an indirect or statistical estimate of the masses is possible. We 
review the general problem of mass measurement in radio pulsar binaries, and 
critically discuss all current estimates of the masses of radio pulsars and their 
companions. We find that significant constraints exist on the masses of twenty- 
one radio pulsars, and on five neutron star companions of radio pulsars. All 
the measurements are consistent with a remarkably narrow underlying gaussian 
mass distribution, m = 1.35±O.O4M . There is no evidence that extensive mass 
accretion (Am > 0.1M©) has occurred in these systems. We also show that the 
observed inclinations of millisecond pulsar binaries are consistent with a random 
distribution, and thus find no evidence for either alignment or counteralignment 
of millisecond pulsar magnetic fields. 
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1. Introduction 

Neutron stars have been the subject of considerable theoretical investigation since long 
before they were discovered as astronomical sources of radio and X-ray emission flBaade a"nd 



Zwicky 1934 , Pppcnhcimcr and Volkoff 1939 , Wheeler 1966| ). Their properties are deter- 



mined by the interplay of all four known fundamental forces — electromagnetism, gravitation, 
and the strong and weak nuclear forces — but neutron stars remain sufficiently simple in their 
internal structure that realistic stellar modeling can be done. Measurements of their masses 
and radii (as well as detailed study of their cooling histories and rotational instabilities) 
provide a unique window on the behavior of matter at densities well above that found in 
atomic nuclei (p nuc ~ 2.8 x 10 14 g cm -3 ). Observations of neutron stars also provide our only 
current probe of general relativity (GR) in the "strong-field" regime, where gravitational 
self-energy contributes significantly to the stellar mass. 

The most precisely measured physical parameter of any pulsar is its spin frequency. The 
frequencies of the fastest observed pulsars (PSR B1937+21 at 641.9 Hz and B1957+20 at 
622.1 Hz) have already been used to set constraints on the nuclear equation of state at high 
densities (e.g., [Friedman et al. 1988|) under the assumption that these pulsars are near their 
maximum (breakup) spin frequency. However, the fastest observed spin frequencies may 
be limited by complex accretion physics rather than fundamental nuclear and gravitational 
physics. A quantity more directly useful for comparison with physical theories is the neutron 
star mass. 

The basis of most neutron star mass estimates is the analysis of binary motion. Soon 
after the discovery of the first binary radio pulsar ( [Hulse and Taylor 1975| ), it became clear 
that the measurement of relativistic orbital effects allowed extremely precise mass estimates. 
Indeed, the measurement uncertainties in several cases now exceed in precision our knowledge 
of Newton's constant G, requiring masses to be quoted in solar units GM Q rather than 
kilograms if full accuracy is to be retained. 

After several recent pulsar surveys, there are now about fifty known binary radio pulsar 
systems, of which five or six are thought to contain two neutron stars. It is thus possible for 
the first time to consider compiling a statistically significant sample of neutron star masses. 
It is our purpose here to provide a general, critical review of all current estimates of stellar 
masses in radio pulsar binaries. The resulting catalog, with a careful, uniform approach to 
measurement and systematic uncertainties, should be of value both to those who wish to 
apply mass measurements to studies of nuclear physics, GR, and stellar evolution, and as 
a guide to the critical observations for observational pulsar astronomers. We begin with a 
discussion of known methods for pulsar mass determination (§@), including a new statistical 
technique for estimating the masses of millisecond pulsars in non-relativistic systems. In §|3] 
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we review all known mass estimates, including new data and analysis where possible. Statis- 
tical analysis of the available pulsar mass measurements is presented in §f|. We summarize 
in §§. 

A second paper will consider mass estimates for neutron stars in X-ray binary systems 
( |Cliakrabarty and Thorsett 19981 , Paper II). A detailed discussion of the implications of 
the combined results of this work and Paper II for studies of supernovae and neutron star 
formation, mass transfer in binary evolution, the nuclear equation of state, and GR will 
occur elsewhere (Paper III). 



2. Methods of mass estimation 



It is a familiar circumstance that estimates of astronomical masses are available only 
for bodies in gravitationally bound binary systems or clusters. Compact stars introduce the 
additional possibility of directly measuring the surface gravitational potential, and hence 
M/R, through the study of redshifted spectral features. Although this technique has been 
used with considerable success in the case of white dwarfs, and attempts have been made to 
fit redshifted X-ray spectra from neutron stars (Paper II), no lines have been identified in 
radio pulsar spectra and other gravitational effects on the observed emission from pulsars are 
sufficiently complex and theory dependent that no useful limits on the neutron star properties 
have yet been possible. In the following, we thus limit ourselves only to the determination 
of stellar masses in binary systems. 



2.1. Pulsar timing 

In any binary pulsar system, five Keplerian parameters can be very precisely measured 
by pulse timing techniques ( [Manchester and Taylor 1977| ): the binary period P&, the pro- 



jection of the pulsar's semimajor axis on the line of sight x = aisini/c (where the binary 
inclination i is the angle between the line of sight and the orbital angular momentum vec- 
tor, defined to lie in the first quadrant), the eccentricity e, and the time and longitude of 
periastron, To and uq. It is frequently more convenient to use the orbital angular frequency 
in place of the orbital period: n = 2tt / P&. These observational parameters are related to the 
pulsar and companion masses, ni\ and m-2, through the mass function 
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where M = mi +m2, the masses are measured in solar units, and we introduce the constant 
Tq = GM Q / c 3 = 4.925 490 947 x l(T 6 s. 

Relativistic corrections to the binary equations of motion are most often parameterized 
in terms of one or more post-Keplerian (PK) parameters (Pamour and Deruelle 19861 , |Taylor| 
and Weisberg 19891 , Pamour and Taylor 1992]) . In GR, the most significant PK parameters 



have familiar interpretations as the advance of periastron of the orbit ill, the combined effect 
of variations in the transverse Doppler shift and gravitational redshift around an elliptical 
orbit 7, the orbital decay due to emission of quadrupole gravitational radiation and the 
"range" and "shape" parameters r and s that characterize the Shapiro time delay of the 
pulsar signal as it propagates through the gravitational field of its companion. In terms 
of measured quantities and the pulsar and companion masses (in solar units), these PK 
parameters are given by ( |'l'aylor 1992| ): 



Co = 3nl(T G M)^ (l-e 2 ) \ (2) 

1 — 4 

7 = en~ s T(|M~ 3 m 2 (jrii + 2m 2 ) , (3) 
192tt s /_ . 73 37 A \ ^ 2 tJ . __i 



P b = —n^l + —e 2 + —e / ')(l-e 2 )-zT3m 1 m 2 M-2, (4) 

r = T m 2 , (5) 
s = xn^T Q 3 M%rri2 l . (6) 



Note, by combining equations ([I]) and (^]), that s = sini for GR.) 

The measurement of the mass function / (equation [l]) together with any two PK param- 
eters (equations @-§D is sufficient in the context of GR to uniquely determine the component 
masses mi and m 2 . With additional assumptions, such as a uniform prior likelihood for or- 
bital orientations with respect to the observer (§ [2.3.4 ), strong statements about the posterior 



distribution of the masses is often possible if even a single PK parameter is measured. 

As an alternative to the PK formalism, which is designed for testing gravitation theory, 
it is sometimes advantageous to fit the timing data to a model that assumes the correctness 
of GR. For example, the DDGR model of Damour, Deruelle, and Taylor (|Taylor 1987| , [Taylor 



and Weisberg 198~9|) describes the pulsar phase as a function of the five Keplerian parameters, 



the companion mass m 2 , and the total mass M. A least-squares fit of the timing data to the 
DDGR model thus gives direct estimates of the uncertainties in m 2 and M, as well as the 
covariances between the mass estimates and estimates of other parameters. 

It is of interest to note that mass measurements from pulsar timing observations depend 
on the unknown relative motion of the solar system and the binary barycenters. Damour 
and Deruelle (1986) have shown that neglecting this velocity is equivalent to changing units 
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of mass and time. In particular, the rest frame mass m and the barycenter frame mass m ssb 
are related by m = Dm ssb , with Doppler factor 

D S '-"-V'/', (7) 

where n is the line of sight unit vector and Vj is the barycentric velocity of the pulsar. 
Although the transverse velocity of the binary system can be estimated from proper motion 
measurements, the radial component n- v 6 is unknown. For a typical velocity of 100 km s _1 , 
the systematic mass error is about ~ 0.03%; small, but in some cases much larger than other 
uncertainties. 



2.2. Masses of companion stars 

While timing measurements of relativistic corrections to the Keplerian orbital equations 
provide the most accurate and theory-independent estimates of neutron star masses, they 
are possible only for close, eccentric binary orbits or when the orbit is observed nearly edge- 
on. In the great majority of observed binaries, the mass function provides the only timing 
information about the component masses, and the pulsar mass can only be determined if 
additional constraints are found on the companion mass and binary inclination through 
other techniques. This section describes several ways to use observations or theoretical 
considerations to limit 777,2; alternate limits on sin i are the topic of § p.3| . 



2.2.1. Optical observations of white dwarf companions 

In recent years, over a dozen companion stars in radio pulsar binaries have been optically 
detected (e.g., van Kcrkwijk 1996 , Lundgrcn, Camilo, and Foster 1996 ). In most cases, the 



companions are white dwarfs; the two main sequence exceptions are discussed in the next 
section. Many of the white dwarf companions are extremely faint (my ~ 26), allowing 
detection only with the Hubble Space Telescope and the largest ground-based telescopes. 

There are a number of ways to determine white dwarf masses (e.g., |Reid 19961 and 
references therein). Given a theoretical relation between the white dwarf mass and radius, the 
measurement of any combination of the mass and radius is sufficient to determine the mass. 
For example, the radius can be estimated directly from estimates of the optical flux, effective 
temperature, and distance. Alternately, the surface gravity, logg, can be found by fitting 
a model atmosphere to the observed spectrum ( [Bergeron, Wesemael, and Fontaine 19911 ), 
and combining the result with a temperature or luminosity estimate. In practice, difficulties 
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arise from several sources. First, the white dwarf companions of millisecond pulsars, with 
typical masses m 2 < O.5M , are usually believed to be helium white dwarfs which were 
insufficiently massive to burn to carbon; the luminosity and temperature evolution of such 
stars has received much less study than for more massive white dwarfs (e.g., |D'Antona and 



Mazzitelli 1990| ), leading to uncertainties in the finite temperature contributions to the mass- 



radius relationship. Hansen and Phinney (1998a,b) have recently calculated cooling curves 
for helium core dwarfs, using their own calculations of H and He opacities at temperatures 
below 6000 K, and applied their results to a number of radio pulsar companions. We discuss 
their results in more detail in ^Q. 

The measurement of surface gravity in cool stars is also potentially problematic. A 
higher helium abundance, as may occur through enhanced convective mixing in cool stars, 
will produce higher pressure and broader lines and hence mimic a higher white dwarf mass 
( [Bergeron, WcscmacL and Fontaine 1991 ); indeed, there is some evidence that surface gravity 



measurements overestimate masses below about 12,000 K ([Reid 1996|) 



2.2.2. Optical observations of main sequence companions 



Two pulsars have been found in binaries with main sequence companions: PSR B 1259— 63, 



with a Be companion flJohnston et al. 1992Q , and PSR J0045— 7319, with a B star companion 
( [Kaspi et al. 1994 ). Although in both cases the optical companion is quite bright and easily 
observed, knowledge of the companion mass m 2 and pulsar mass function / is of little use in 
limiting m\ when m\ <C m 2 . If the mass function of the companion can also be measured, 
then the mass ratio can be determined. This has been done in the case of J0045— 7319 

OCT- 



2.2.3. The P&^m 2 relation 

The binary millisecond pulsars are believed to be spun up to high spin frequencies 
through mass transfer from a companion star. In many cases, these pulsars are observed in 
wide, low-eccentricity binary systems with white dwarf secondaries. These characteristics 
indicate that the secondary must have passed through a red giant phase after the forma- 
tion of the neutron star primary, during which tidal torques circularized the orbit and the 
giant probably filled its Roche lobe, causing the mass transfer which spun up the pulsar to 
millisecond periods. At the end of mass transfer, the envelope of the giant is exhausted or 
ejected, leaving the degenerate core as a white dwarf secondary. There is a close relation 
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between the core mass and the radius of low mass giants ( [Refsdal and Wcigcrt 1971 , Web 



bink, Rappaport, and Savonije 1983| , |Joss, Rappaport, and Lewis 1987| , |Rappaport et al. 



1995| , |Rappaport and Joss 1997| ). Combined with the assumption that the giant filled its 
Roche lobe during the mass transfer, this yields a relation between the binary period at the 
end of mass transfer and the remnant white dwarf mass ( [Rappaport et al. 1995| ) . 

Following Rappaport et al. (1995) and references therein, the relation between the ef- 
fective radius of the Roche lobe Rl and the binary separation a can be written 



R L w 0.46 a 1 + — 



\ -1/3 

mi \ 



m 9, 



where m g is the total mass of the giant (core and envelope). Near the end of mass transfer, 
the envelope mass can be neglected, so that m g m 2 (where m 2 is the final white dwarf 
mass), and the giant radius R g ~ Rl- Using Kepler's third law, 

P b = 0.374 R z J 2 mz 1/2 days, (9) 

where R g and m 2 are in solar units. We note that equation (|9]) is independent of the pulsar 
mass rrii, and relates the orbital period at the end of mass transfer to the final white dwarf 
mass and the giant radius R g . The radius depends, in general, on the composition and 
history of the giant as well as on m 2 ; the utility of equation (|]) derives from the relatively 
narrow distribution of R g for a given m 2 . 

Rappaport et al. (1995) have examined stellar models with a wide range of chemical 
compositions, varying between Population I and Population II values, as well as a variety of 
envelope masses and convective mixing lengths. Over a wide range of core masses greater 
than O.15M , they find the data are well described by the equation 

*» = TT4^ + °' 8 (10) 

where _R — 495Oi? is the best fit to the stellar models. In all models studied by Rappaport 
et al. , equation ( [Of) was correct within a factor of 1.8. Over the limited range m 2 < O.25M , 
Rappaport and Joss (1997) found that the alternate expression 

log Rg = 0.031 + 1.718 m 2 + 8.04 ml (11) 

could be used, with the smaller uncertainty of a factor ~ 1.3. We regard these to be 
approximately 95% confidence regions for R g . 

It is important to determine the range of applicability of equations ( Jl~0| ) and (p4J). At 
orbital periods less than a few days, heating of the surface of the companion by X-rays 
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produced during accretion may cause significant bloating and modification of the core mass- 
radius relations ( |Podsiadlowski 1991| , [Rappaport et al. 1995| ). Following Rappaport et al., 



we trust the modeling only for binary orbits longer than three days. Further, the Pb-m 2 
relation cannot be applied to systems like PSR J2145— 0750 or J1022+10 — both 16 ms pulsars 
that were most likely recycled in common envelopes that occurred when their companions 
overflowed their Roche surfaces while on the asymptotic giant branch, leaving carbon-oxygen 
rather than helium white dwarfs (|van den Heuvel 1994|) — or to pulsars like B0820+02, a slow 



(0.864 s), high-field (3 x 10 11 G) pulsar without evidence of significant recycling. To avoid 
both classes of pulsar without introducing biases by cutting on apparent companion mass, 
we apply the P\ ) -m2 relation only to pulsars with P < 10 ms, for which the assumption of 
an extended period of mass transfer during a low-mass X-ray binary phase seems secure. 

Using equations (P|)-(|Tl"|), the observed orbital period in a millisecond pulsar binary can 
be used to limit the range of rri2- When combined with the mass function (equation HJ) and 
the restriction sini < 1, an upper limit can be placed on the pulsar mass. The resulting 



limits are discussed below, in §3.2j . If a random distribution of orbits is assumed for the set 
of observed millisecond pulsars, then statistical arguments lead to limits on the masses of 
milliseconds pulsars as a class. These arguments are described in §01. 



2.3. Alternative methods of estimating the orbital inclination 

If the companion mass can be estimated and bounds can be found on the orbital incli- 
nation, then the mass of the pulsar can be found using the Keplerian mass function. In this 
section, we discuss measurements that may allow direct estimates of sini. 



2.3.1. Polarization measurements 

In the standard model of millisecond pulsar formation, the pulsar is spun-up by mass 
transfer from a companion star. After spin-up, the pulsar spin axis will be aligned with the 
orbital angular momentum, so a measurement of the angle ( between the pulsar spin axis 
and the line of sight also determines the orbital inclination: i = (, or i = tt — ( for ( in the 
second quadrant. (Although this is also true at an early stage of the formation of double 
neutron star binaries, an asymmetry in the second supernova explosion may leave the spin 
and orbit misaligned.) 

In the rotating vector model of pulsar emission, the radio signal is elliptically polarized 
with the major axis of the polarization ellipse aligned with the plane of curvature of the 
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dipolar magnetic field lines. The position angle of the observed linear polarization ifj can be 
expressed as a function of the rotational phase of the pulsar: 

sin a sin 6 , , 

tan V# - V>o = -t-t t, (12) 

sin (, cos a — cos (j sin a cos 

where a is the angle between the pulsar spin axis and the magnetic pole. In practice, the 
difference ( — a can often be estimated quite well from equation flT2|), but accurate estimates 
of ( or a are possible only if polarized emission can be detected over a broad range of pulse 
phase <p. 

Application of equation fll2j) to millisecond pulsars, in which the light cylinder bounds 
the magnetosphere at only a few stellar radii, may be complicated by multipolar field ge- 
ometries, aberration, and other deviations from the simple rotating vector model. Indeed, 
attempts to fit equation ( |I2] ) to millisecond pulsar data have met with few unqualified suc- 
cesses ( [Thorsett and Stinebring 1990| , [Navarro et al. 19971 ). 



2.3.2. Interstellar Scintillation 

Observed pulsar signal strengths vary in both frequency and time because of scintilla- 
tion in the interstellar medium. A phase- changing screen along the line of sight produces 
an interference pattern across which the pulsar moves; the characteristic scintillation decor- 
relation timescale ri ss is inversely proportional to the transverse component of the pulsar 
velocity. Observations of scintillation rates have been widely used to estimate the proper 
motions of isolated pulsars (e.g., |Cordes 1986|) . In binary systems, the transverse velocity 



is modulated by the orbital motion, with a small amplitude of modulation for orbits viewed 
face-on, and a large amplitude for orbits viewed on edge. Scintillation measurements can 
thus provide an alternate path to estimating the inclination angle sini ( |Lyne 1984|) . 

The real physical situation is more complex. The measured scintillation velocity v iss of 
a pulsar depends not only on its proper motion v pm and orbital motion v orb (t), but also on 
the Earth's motion v^(t), the mean motion of the scattering medium Vj sm and the ratio of 
the effective distance to the scattering screen to the pulsar distance, /: 

Vi ss (t) = (1 - /) (V rb(*) + Vpm) + /v® (t) - V ism . (13) 

If the proper motion is known, from timing measurements or interferometry, then the annual 
modulation of and the orbital modulation of v or b can in principle be used to find not 
only the orbital inclination sini but also /, the pulsar distance d, and the position angle on 
the sky of the orbital ascending node Q. 
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Only a few attempts to determine orbital inclinations from scintillation observations 
have been published. In no case was v pm known, nor has the annual modulation due to 

been measured (although the latter has been observed in the isolated millisecond pulsar 
B1937+21 by Ryba (1991)). Lyne (1984) found that the inclination of PSR B0655+64 is 
either 62° or 84° (the discrete ambiguity could be broken through observations at another 
time of year, using variations in vm). Jones and Lyne (1988) later expressed reservations 
about this measurement. In the most convincing success of the scintillation technique, Dewey 
et al. (1988) found a limit on the the inclination of PSR B1855+09: sin 2 > 0.94. Shapiro 
time delay measurements have since shown that sini rs 0.9993 (§ 3.2.7[) . 



The problem of estimating scintillation parameters from data is considered by Cordes 
(1986). The dominant error source in estimating scintillation parameters is the finite number 
iV of scintillation features sampled: a v /v- lss « 0.6/ y/N. Scintillation intensity fluctuations 
are exponentially distributed; for observing bandwidth B and time T, we have roughly 
N ~ 10~ 2 i?T/(Ai/ iss Ti SS ). Considering a case relevant to many recently discovered pulsar 
binaries, a pulsar at a dispersion measure DM= 20 will have typical scintillation parameters 
at 430 MHz r iss ~ 10 minutes and Az/ iss ~ 100 kHz QRickett 1988|) . With B = 10 MHz 
and T = 1 hr, <r v /v iss ~ 25%, so a study of orbital dependence of scintillation parameters 
requires a substantial amount of observing time. Furthermore, to use a bandwidth of 10 MHz 
requires a spectrometer with <; 100 frequency channels. Such observations will be more easily 
done with the new generation of flexible, all-digital pulsar data recorders ( |Shrauner et al.\ 
1996| , |Jenet et al. 1997|) . It is possible that long term variability in / or v iss will limit the 



use of the annual variation in vm to estimate /. 



2. 3. 3. Secular variation of x = a\ sin i 



Proper motion of the binary system across the sky leads to a secular change in the 
projected semimajor axis x. If Q is the position angle of the ascending node and \i a and \i§ 
are the components of the proper motion \l in right ascension and declination, then 



x 

— = cot % 

X 



-fjL a sin Q + us cos Q) . 



(14) 



The angle Q is generally unknown (though it is accessible to scintillation measurements, 
§ 2.3.2j ), but eqn. can be rewritten as the limit 



tanz < 



x 
x 



(15) 



A similar expression can be found relating proper motion to a purely geometric advance of 
periastron Co — in principle, measurement of both x and Co would determine i — but the effect 
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is too small to have yet been seen in any pulsar binary. 



2.3.4- Random distribution of orbital inclinations 

When no other information is available on the orbital inclination, it is sometimes useful 
in statistical calculations to assume that binary orbits are randomly oriented with respect to 
the line of sight. The differential distribution of inclinations is then proportional to sinz (i.e., 
the most likely orbital viewing angle is edge-on). Values of cosi should then be uniformly 
distributed between and 1. 

In binary systems where the directions of the pulsar spin axis and the orbital angular 
momentum vector are expected to be correlated, such as millisecond pulsar systems in which 
significant mass transfer has occurred since the last supernova, non-random orientations of 
the pulsar magnetic field axis with repect to the spin axis may lead to a pulsar discovery bias 
that skews the observed orbital orientation distribution. Some models of pulsar evolution 
predict magnetic field alignment, counteralignment, or both (e.g., |Ruderman 199 1| ) . 



Backer (1998) has claimed, based on the distribution of observed mass functions in 
21 pulsar-white-dwarf binaries, that millisecond pulsars have a non-random distribution of 
observed binary inclinations. He suggests that an apparent preference for high inclination 
orbits may arise because the pulsars' magnetic fields are preferentially oriented perpendic- 
ular to the spin axis. However, his analysis makes several assumptions and approximations 
that prove unwarranted. In particular, he compares the observed distribution of minimum 
companion masses m^mm with the predicted distribution given the observed mass functions, 
random inclinations and a single fixed value for the true companion masses m2 (i.e., a delta 
function distribution in m 2 ), finding poor agreement. However, the predicted distribution 
for m 2im in actually depends quite sensitively upon the assumed distribution for m 2 ; the delta 
function distribution is an unreasonable assumption. Indeed, the -Pf,-m 2 relation suggests 
that the observed m 2 values should vary by a factor of three. Knowing the precise limits 
and distribution of these masses is crucial for analyzing the distribution of m 2im i n . 

We repeated Backer's analysis, with different assumptions about the distribution of m 2 
values. If we assume the true values m 2 fall uniformly in the range 0.15-0.35 M , then the 
predicted and observed distributions for m 2 min agree at an 81% confidence level according to 



a Kolmogorov-Smirnov (KS) test (e.g., [Eadie et al. 19711) . However, assuming a 0.1-0.4 M ( 



' G 



range reduces the agreement to the 33% confidence level. The discrepancy is primarily due to 
a deficit of observed systems with small values of m 2 min , as noted by Backer. Unfortunately, 
our knowledge of the true underlying distribution of m 2 in these binaries is sufficiently 
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uncertain that no conclusion regarding the distribution of inclinations is possible from this 
line of analysis. However, we argue below (§ffl) that for systems with orbital periods longer 
than three days, where the -P{,-m 2 relation can be applied to estimate the m 2 distribution, 
the observed mass functions are consistent with a random distribution of inclinations. 



3. Radio pulsar masses 

The parameters of 47 radio pulsar binary systems are listed in Table [1]. At least five 
of these pulsars have neutron star companions, so a total of 52 neutron stars are known in 
radio pulsar binaries. Useful mass constraints can be placed on half of these stars. 



3.1. Double neutron star binaries 

3.1.1. PSR J 15 18+ 4904 

PSR J1518+4904 is in a moderately relativistic binary system; one PK parameter has 
been measured, the relativistic advance of periastron Co = 0.0111(2)°yr~ 1 , yielding a total 



system mass of 2.62(7)M ([Nice, Sayer, and Taylor 1996|) . Given the mass function / = 
O.115988M , the lower limit on the companion mass (using sin 2 < 1) is m 2 > (/M 2 ) 1 / 3 = 
O.93M , and the lower limit on the inclination (given mi > 0) is sinz > (f/M) 1 / 3 = 0.35 
(or i > 20°). Using a uniform prior distribution in sini over the interval 0.35 < sin 2 < 1, 
we find the central 68% confidence intervals mi = 1.56l[J;44 an d m 2 = l-OS^-ii an d central 
95% confidence intervals mi = l.SG^Jo and m 2 = l-05ljjj;?4. 

The identification of the companion as a neutron star is compelling, although it is not 
possible to completely rule out a low-mass black hole. Optical observations (van Kerkwijk, 
quoted in Sayer 1996) show no source at the pulsar position to ~ 24.5, excluding a main 
sequence companion, and a white dwarf in an eccentric orbit is not expected on evolutionary 
grounds. 



3.1.2. PSR B1534+12 



PSR B1534+12 is in a highly relativistic binary system, presumably with a second 
neutron star companion. All five of the PK parameters defined by equations (f2|)-(|6]) have 
been measured, three with better than 1% precision ( |Stairs et at. 1998"D - Using the DDGR 
timing model, the total system mass is found to be M = 2.67838(8)M Q , while the individual 
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component masses mi and m 2 are both 1.339(3)M . (The quoted errors are 68% confidence 
regions; the 95% confidence regions are about twice as large.) It is remarkable that the pulsar 
and companion masses agree to better than 0.4%; the assumption that the companion is a 
second neutron star seems secure. The uncertainty on the individual masses is expected to 
decrease significantly in the next year, when observations are again possible with the Arecibo 
telescope. 



3.1.3. PSRB1913+16 

PSR B1913+16 was the first binary pulsar discovered and, with observations stretching 
over two decades, it remains one of the best studied. PSR B1913 is a highly relativistic 
system: the PK parameters ui, 7, and Pb have all been measured precisely. A fit to the 
DDGR timing model yields the total system mass M = 2.82843(2) and component masses 
mi = 1.4411(7) and m 2 = 1.3874(7)M ( [laylor 1992]) . (The quoted errors are 68% confi- 



dence regions; the 95% confidence regions are about twice as large.) The uncertainties are 
approaching the level where they will be dominated by kinematic effects (equation |7|). 



3.1.4. PSR B2127+11C 

PSR B2127+11C is a relativistic binary in the globular cluster M15. Precise mea- 
surements have been made of the PK parameters Co and 7, resulting in mass estimates for 
the pulsar and companion of 1.349(40) and 1.363(4O)M , respectively, and a total system 
mass M = 2.7121 (6)M Q ( peich and Kulkarni 1996] ). (The quoted errors are 68% confidence 
regions; the 95% confidence regions are approximately twice as large.) 



3.1.5. PSR B 2303 +46 

PSR B2303+46 is a pulsar in a moderately relativistic orbit with, most likely, a neu- 
tron star companion. One PK parameter, the relativistic advance of periastron, has been 
measured. A timing analysis was published by Thorsett et al. (1993), and updated by Ar- 
zoumanian (1995). Using data from early 1985 to late 1994, we find an improved value 
Co = 0.01019(13)°yr _1 , yielding a total system mass M = 2.64 ± 0.05M Q and the constraints 
mi < 1.44M and m 2 > 1.2OM . Using a uniform prior distribution in sini, we find the cen- 
tral 68% confidence intervals mi = 1.30lo; 4 g an d m 2 = ^-^-o.il an d central 95% confidence 
intervals m 1 = 1.30lS$j and m 2 = 1.34±J;?|. 
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3.2. Neutron star/white dwarf binaries 

3.2.1. PSR J0437-4715 

PSR J0437— 4715 is the brightest and closest known millisecond pulsar, and therefore one 
of the best studied. It has a mass function / = 1.243 x 10~ 3 M Q , and a white dwarf companion 
that has been detected optically flJohnston et al. 1993| , [Bailyn 1993| , [Bell, Bailes, and Bessell 



1993| , |Danziger, Baade, and Delia Valle 1993|) ; and both thermal and nonthermal X-ray 



emission from the neutron star have been observed (Pecker and Triimper 1993] , |Halpern 



[Martin, and Marshall 1996[ ) . The distance is known, d = 178 ± 26 pc, from the effects of 



parallax on the pulsar timing signal flSandhu et al. 1997] ). Using the optical data of Danziger 
et al., Hansen and Phinney (1998b) found an effective temperature of the companion T c g = 
4600 ± 200 K. Combining the T e g and distance estimates with their cooling curves, they find 
consistent companion models for all masses 0.15 < 7772 < O.375M . This range encompasses 
the mass limits derived from the -P&-7772 relation (§ 2.2.3 ), 0.16 < m 2 < 0.23M Q . 



The binary orbit is nearly circular, and despite the very high timing precision achieved, 
no PK timing parameters have been measured. However, the proximity of the pulsar leads 
to a high proper motion (/i = 141 mas/yr), changing the projected orbital size ( §2.3.3] ) at a 



rate x/x = 2.43(12) x 10 14 flSandhu et al. 19TJ7 ). The implied limit on the inclination angle 
is i < 43°, or sin 2 < 0.682. 

Attempts to measure the system geometry using the polarization of the radio beam 
(§ |2.3.1 ) have proven difficult, because of the very complex emission pattern. Well-calibrated 



intensity and polarization data have been reported by Manchester and Johnston (1995), who 
have modeled the sweep of the linear polarization position angle across the profile in terms 
of the rotating vector model. They find that an impact parameter of the line of sight on the 
magnetic pole (3 = — 5° and an angle between the line of sight and the spin axis ( = 140° 
produced reasonable agreement with the data, but there are strong systematic deviations 
from the simple rotating vector model. Using the same data, Gil and Krawczyk (1997) model 
the multicomponent profile and find a comparable impact parameter, (3 = —4°, but a very 
different £ = 16°. They show that this geometry also explains the observed polarization 
sweep over much of the pulse period. With the assumption that the pulsar spin axis is 
aligned with the orbital axis, the implied orbital inclination is sin 2 = 0.64 (Manchester and 
Johnston model) or sini = 0.28 (Gil and Krawczyk model). Either is consistent with the 
timing data. 

The timing limit on sinz, together with the mass function and the upper limit on 
m2 from the core mass-orbital period relation, gives an upper limit to the pulsar mass of 
mi < 1.51M , while the optical upper limit on mi gives only the much weaker limit m\ < 
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3.29M . The Manchester and Johnston inclination together with the Pb-m 2 mass range gives 
0.77 < m 1 < 1.37M , while the Gil and Krawczyk inclination gives 0.11 < m x < O.23M , a 
physically implausible result. In any case, the complexity of the pulse shape modeling leads 
us to adopt the weaker one-sided limit mi < 1.51M . 

Improvements of the parallax measurement and the optical photometry and spectroscopy 
are needed to further constrain the companion mass. An independent mass determination 
would test the core mass-orbital period relation as well as improve the limit on m\. Mea- 
surement of the position angle of the ascending node (perhaps through scintillation studies), 
combined with the measurement of x/x, would give the inclination of the orbit. 

3.2.2. PSR J 1012+5307 

PSR J1012+5307 has a hot, bright white dwarf companion that has been extensively 
studied, van Kerkwijk, Bergeron, and Kulkarni (1996) have used the models of Bergeron, 
Wesemael, and Fontaine (1991) to determine the effective temperature T e g = 8550 ± 25 K 
and surface gravity log g = 6.75 ±0.07 of the companion. The latter value is in disagreement 
with an unpublished value of Callanan and Koester, log g = 6.4 ± 0.2, quoted by Hansen 
and Phinney (1998b). Using the optical observations and their cooling models for low mass 
helium white dwarfs, Hansen and Phinney find a companion mass 0.165 < m 2 < O.215M 
for the van Kerkwijk et al. gravity measurement and 0.13 < m 2 < 0.18 for the Callanan and 
Koester value. 

The radial velocity of the companion has been measured, making J1012+5307 a double- 
line spectroscopic pulsar binary, van Kerkwijk et al. have found mi/m 2 = 9.5±0.5 (van Kerk- 
wijk, private communication; their earlier published value mi/m 2 = 13.3 ±0.7 was corrupted 
by a calibration problem). Depending on the gravity measurement, the resulting (la) pul- 
sar mass limit is 1.5 < mi < 2.2M or 1.2 < mi < 1.8M . Clearly resolution of the 
discrepant gravity measurements must be a high priority, the radial velocity measurements 
now contribute negligibly to the total error on m 2 . For now, we adopt Hansen and Phinney's 
conservative conclusion that 1.2 < mi < 2.2M , and we regard the error range as roughly a 
68% confidence interval. 

3.2.3. PSR J 1045-4509 

PSR J1045-4509 is in a 4.08 day orbit with mass function 1.765 x 1O _3 M . The P b -m 2 
relation gives an upper limit to the companion mass of O.168M , leading to a limit on the 
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pulsar mass of mi < 1.48M . 

3.2.4. PSR J 1713+0747 

PSR J1713+0747 is a bright pulsar in a nearly circular, 68 day orbit with a white dwarf 
companion. At the timing precision reached (about 500 ns for 1.4 GHz observations), Camilo 
(1995) found that the Shapiro PK parameters r and s were required to adequately model the 
data, but strong covariances between the range parameter r and other orbital parameters 
(especially x) prevented him from setting interesting limits on the pulsar or companion mass. 

The P\,-vn,2 relation predicts 0.26 < m 2 < 0.35M Q . Companion mass estimates are also 
possible by combining unpublished optical observations of the white dwarf by Lundgren et al. 
with the parallax from pulsar timing flUamilo, Foster, and Wolszczan 1994| ). The resulting 
limits are 0.15 < m 2 < O.31M if the white dwarf has a thick H envelope, and m% < O.27M 
if it has a thin H envelope ( |Hansen and Phinney 1998b| ). 

If the observed mass function / = 7.896 x 1O _3 M is combined with the Pfo-m 2 limit on 
m 2 and the restriction sini < 1, an upper limit to the pulsar mass is found: m% < 1.94M . 
The optical observations yield a tighter limit, m\ < 1.63M , as was noted by Hansen and 
Phinney. More interesting limits can be obtained by combining the Shapiro measurements 
with the limits on m 2 (recalling that m 2 = r). For example, if m 2 = O.299M , the timing data 
of Camilo is consistent with only the restricted range sin i = 0.963(4). By letting m 2 vary over 
the range allowed by the -P&-m 2 relation, we find the allowed pulsar mass m 1 = 1.45 ± 0.31, 
where the uncertainties are dominated by the systematic uncertainties in the -P&-m 2 relation. 
If the companion mass is limited to the intersection of the regions allowed by Pf,-m 2 and the 
optical measurements, 0.26 < m 2 < O.31M , then we find mi = 1.34 ± O.2OM . 

3.2.5. PSR B 1802- 07 

PSR B1802— 07 is in the globular cluster NGC 6539. Its companion is most likely a 
white dwarf; the system's large eccentricity can be understood as the result of gravitational 
perturbations of the system by close stellar encounters in the dense cluster. 

The relativistic advance of periastron Cj was measured by Thorsett et al. (1993), and 
an improved measurement was published by Arzoumanian (1995). Using the same analysis 
techniques, with data extending through October 1997, we find a slightly improved value 
Co = 0.0578(16), implying a total system mass M = 1.62(7)M . The lower bound on the 
companion mass from the requirement that sinz < 1 is (/M 2 ) 1 / 3 = O.29M , and the lower 
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bound on the inclination from the requirement mj > is sini > (//M) 1 / 3 = 0.18 (or % > 10°). 
Using a uniform prior for sini in the range 0.18 < sin 2 < 1, we find a 68% confidence bound 
on the pulsar mass m 1 = 1.26io;$fM Q and a 95% confidence bound mi = 1.2611q'^Mq. 



3.2.6. PSR J1804-2718 

PSR J1804— 2718 is in an 11.1 day orbit. The P^-m^ relation limits its companion mass 
to 0.185 < rri2 < O.253M ; the resulting upper limit on the pulsar mass is mi < 1.73M Q . 



3.2.7. PSR B '1855 +09 

PSR B1855+09 is a 5.4 ms pulsar in a 12.3 day circular orbit (e = 2 x 10~ 5 ) with a white 
dwarf companion. Because it is a bright pulsar for which high precision timing measurements 
can be made and because the orbital inclination is high, measurements have been made of 
the PK parameters r and s. 

We have reanalyzed all of the available data, extending from January 1986 to January 
1994. (For the data collected after mid-1989, taken with the Princeton Mark III timing sys- 
tem ( [Stinebring et al. 1992]) , a new algorithm was used for identifying times and frequencies 
where the signal strength was enhanced by interstellar scintillation and weighting this data 
in subsequent analysis.) The details and complete timing solution will be published else- 
where. The Shapiro parameters are r = Gmj/c 3 = 0.248(1 1)M Q and s = sini = 0.9993(2). 
The resulting limits on the pulsar mass are mi = 1.41(1O)M (68% confidence; the 95% 
confidence region will be about twice as large). The new measurement is in good agreement 
with previous values: m x = 1.27±o;i 3 M QRyba and Taylor 1991|) and m\ = 1.50io]i4 QKaspi 



Taylor, and Ryba 1994Q . Q 



The P;,-m2 relation predicts 0.19 < ui2 < O.26M , in good agreement with the timing 
measurement. The corresponding upper limit to the pulsar mass is m x < 1.51M Q . As noted 
above, scintillation measurements have also been used to limits sini > 0.94 ( |lJewey et al 
19881) . 



3 Note that both Ryba and Taylor (1991) and Kaspi et al. (1994) incorrectly (over) estimate the error on 
mi by calculating the joint 68% confidence bound on r = T^rn-i and s — sini and then interpreting this as 
a 68% confidence bound on each of rri2 and sini. 
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3.2.8. PSR J2019+2425 

PSR J2019+2425 is in a 76.5 day orbit. The -P&-m 2 relation gives a limit on the 
companion mass 0.264 < m 2 < O.354M ; the resulting upper limit on the pulsar mass 
is mi < 1.68M . 



3.3. Neutron star/main sequence binaries 

3.3.1. PSR J 0045- 7 3 19 

The pulsar J0045— 7319 is the only known pulsar in the Small Magellanic Cloud. It is in 
a binary orbit, with mass function / = 2.17M ( [Kaspi et al. 1994]) . The companion has been 



identified as a Bl V star. The radial velocity of the companion has been measured, giving a 
mass ratio q = m 2 /mi = 6.3 ± 1.2 ( Pell et al. 1995|) We have compared the observed optical 



luminosity L = 1.2 x 1O 4 L and temperature T e g = 2.4(1) x 10 4 K to the grids of stellar 
models calculated by Schaller et al. (1992) for the low metallicity {Z = 0.001) appropriate 
for the SMC, and estimate the companion mass to be 10 ± 1M . The pulsar mass is then 
mi = m 2 /q = 1.58 ± O.34M (68% confidence), where the uncertainty is dominated by the 
uncertainty in the amplitude of the companion's radial velocity curve. 



4. Discussion 

For a dozen neutron stars, useful mass constraints are available with no assumptions 
beyond the applicability of the general relativistic equations of orbital motion to binary 
pulsar systems. Ten of these stars are members of double neutron star binaries. With 
the possible exception of PSR B2127+11C, in the globular cluster M15, the pulsar in each 
system is believed to have undergone a short period of mass accretion during a high-mass 
X-ray binary phase (Am ~ 1O _3 M , Taam and van den Heuvel 1986). The companion stars 
have not undergone accretion; their masses most directly preserve information about the 
initial mass function of neutron stars. 

Only two "millisecond" pulsars, the end products of extended mass transfer in low-mass 
X-ray binaries, have interesting mass estimates based on GR alone: PSRs B1802— 07 and 
B1855+09. Because such pulsars must accrete ~ O.1M to reach millisecond periods (Taam 
and van den Heuvel 1986|) , and much more (~ O.7M ) in some field decay models (e.g., 
[van den Heuvel and Bitzaraki 1995| ), obtaining additional mass measurements of millisecond 
pulsars is of particular interest in testing evolutionary models and in locating the maximum 
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neutron star mass. 



As noted in § [2.2.3 , the p,-m 2 relation can be used to estimate the companion mass in 



recycled binary systems with circular orbits and orbital periods P& ^ 3 d. There are now 
thirteen such millisecond (P sp in < 10 ms) pulsars known, excluding those in globular clusters 
(where gravitational interactions may have significantly perturbed the orbital parameters 
since spin-up). In each case, the measured mass function and the inferred companion mass, 
together with the requirement that sini < 1, then yields an upper limit on the mass of the 
pulsar itself. A number of systems in which this upper limit is particularly constraining have 



been mentioned in 01372 . 



Additional constraints on the neutron star mass in these systems can be derived using 
statistical arguments, given a prior assumption about the distribution of binary inclinations. 
The simplest such assumption is that the binaries are randomly oriented on the sky, though 
biases toward high or low inclinations are possible in some models (§ |2.3.4| ). However, as 
discussed below, we believe there is currently no evidence for such a bias, so for the remainder 
of this discussion we assume random orbital orientations. 

For an individual system, we are interested in the probability distributionQp(mi; /, P&) 
for the neutron star mass mi given the measured mass function / and binary period P&. We 
can neglect the measurement uncertainty in / and P&. Then, the probability distribution for 
mi can be written schematically as 

p(mi;f,P b )= / d(cosi) c?m 2 p(w 2 ; P&)p(cosi)p(mi|m 2 , cosi; /), (16) 

JO JnT>Z,min(Pb) 

where mi — /-^(masinz) 3 / 2 - m 2 is restricted to positive values. We have evaluated this 
numerically for each system, assuming that p(cosi) is uniform between zero and unity and 
that p(m 2 ; Pb) is uniformly distributed within the appropriate factor (see § |2.2.3| ) of the m 2 



implied by equations (9)— (1 1). Not surprisingly, the width of the distribution p(mi;/, P&) 
is dominated by the range of allowed cosi rather than the uncertainty in m 2 for a given 
Pb- For each of the 13 binaries, we have plotted the cumulative distribution CDF(mi) = 
Jo™ 1 P{ m 'i) dm! x in Figure |l|. 

The median and 68% and 95% confidence regions for each pulsar mass is given in Table ^. 
Although several of the pulsars have, under the assumptions made, most likely masses well 
above 2M Q , some such results are expected even if all the masses are quite low. In fact, in 



4 We adopt the notation that p(x; A) is the (marginal) probability density for the random variable x, 
where x depends upon the parameter A. Also, p{x\y\ A) is the conditional probability density for the 
random variable x for a given value of the random variable y and parameter A. 
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only one case of the 13 pulsars does 1.35M lie outside the 95% central confidence region 
(J1045— 4509), and in 6 cases of 13 is 1.35M excluded at 68% confidence, consistent with 
chance. 

It is interesting to ask whether a single, simple distribution of neutron star masses 
is consistent with all of our observational constraints. We considered two models for this 
question: a Gaussian distribution of masses with mean m and standard deviation a, and a 
uniform distribution of masses between mi and m u (cf. Finn 1994). A maximum likelihood 
analysis was used to estimate the parameters m, a, mi, and m u (assuming a uniform prior 
distribution for all four parameters). The resulting 68% and 95% joint confidence limits on 
m and a are shown in Figure 0, and on m\ and m^ in Figure [| In each model, the 
distribution of neutron star masses is remarkably narrow: the maximum likelihood solutions 
are m = 1.35M© and a = O.O4M , and mi = 1.26M and m u = 1.45M . 

Of course, any model (even a poor one) will yield maximum likelihood parameters 
for a given data set. However, it is obvious by inspection that both the Gaussian and 
uniform distributions for the neutron star mass are good fits to the extremely narrow observed 
range of neutron star masses in the double neutron star binaries. While it is difficult to 
quantify the goodness-of-fit for the entire data set, because of the diverse assumptions made 
in the various mass estimates and the sometimes highly non-gaussian error estimates, we can 
easily test some neutron star subsamples against the maximum likelihood gaussian model 
mi = 1.35 ± O.O4M . For the thirteen neutron-star-white-dwarf binaries, we used a Monte 
Carlo technique to evaluate the fit quality. For each binary (with its measured and /), 
we simulated a large number of Monte Carlo trials where the neutron star mass m\ was 
drawn from the maximum likelihood model, m<i was drawn from the appropriate uniform 
distribution implied by Pf>, and cosi was drawn from a uniform distribution. The Monte 
Carlo trials were then used to construct the probability distribution for the mass function, 
and this distribution was used to compute the cumulative probability for the measured mass 
function, p(f < /). If the model and the associated assumptions are correct, then the 
cumulative probabilities for the 13 measured mass functions should be consistent with a 
uniform distribution betwen zero and unity (cf. a classical V/V max test, Schmidt 1968). A 
KS test of the distribution shows consistency with a uniform distribution at the improbably 
good 99% level (Figure |). 

For the ten stars for which gaussian error estimates a e are available (both stars in 
the relativistic binaries B1534+12, B1913+16, and B2127+11C, as well as J1012+5307, 
J1713+0747, B1855+09, and J0045-7319), we can calculate a \ 2 statistic, £(m-m) 2 / '(a 2 + 
a 2 ) = 7.5, consistent with expectations for a chi-square distribution with 10 — 2 degrees of 
freedom. 
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We conclude, therefore, that at least in the radio pulsar systems, there is no evidence 
for neutron star masses above about 1.45M . Indeed, the data appear very well modeled by 
very narrow distributions centered around 1.35M . 

5. Summary 

There are now 26 neutron stars in binary radio pulsar systems for which useful mass 
constraints can be derived. Of these, about half are neutron star-white dwarf binaries in 
which the mass determination depends on the validity of the P^-wii relation and the isotropy 
of the binary orbits with respect to the line-of-sight, as discussed in the previous section. 
All other mass constraints are listed in Table ^| and shown in Figure ||. 

Although we defer a full discussion of the underlying neutron star mass distribution 
and it implications for neutron star formation and evolution until after analysis of the X-ray 
binary systems (Paper II), we note here a few points of particular interest about the radio 
pulsar binaries. Figure |5] is striking primarily for the very small variations in the masses of 
well measured stars. This has, of course, been noted before (e.g., Thorsett et al. 1993), but 
it remains surprising that no new mass measurements differ greatly from 1.4M Q . In the five 
double neutron star binaries for which a relativistic periastron advance yields an accurate 
total (and hence average) mass, the average neutron star masses vary by less than 7%. 

The most surprising implication of the current results is that there is little evidence for 
mass tranfer of O.1M or more in the millisecond pulsar systems. It is important, then, to 
reiterate the assumptions upon which this conclusion rests: (1) the -Pft-m 2 is correct, at least 
within the (modest) claimed precision, and (2) binary orbits are randomly oriented with 
respect to the line of sight. There are good prospects for testing both assumptions. 

The reliability of the -Pfe-m 2 relation depends principally upon the core-mass-radius 
relation for red giants (eqns. [H] and O). The latter relation can be tested observationally by 
careful study of nearby red giants. Precise measurements of bolometric flux and angular size 
(through optical/infrared photometry and interferometry; see, for example, Dyck et al. 1996, 
Perrin et al. 1998) together with accurate distance measurements (e.g., using Hipparcos) 
can be used to probe the relationship between luminosity and radius. Since the core-mass- 
luminosity relation for red giants should be inherently more precise than the core-mass-radius 
relation (Rappaport et al. 1995 and references therein), such observations would provide an 
effective test of the core-mass-radius (and hence Pft-m 2 ) relation. Further, we hope that 
the power of the Pfe-m 2 relation as a statistical tool will encourage more detailed theoretical 
investigations of, in particular, the extension of the relation to orbital periods below about 
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three days, where X-ray heating and bloating of the companion star become important. 



Although any technique that limits sin 2 (§ |2.3| ; e.g., polarization or scintillation studies) 
can be used to test the hypothesis that orbits are randomly inclined, the most precise mea- 
surements will come from pulsar timing and Shapiro time delay measurements. In Table |2] 
we list the mean predicted amplitude of the Shapiro delay signal At$ = 2m 2 T Q log(l — sini) 
for each of the thirteen pulsars discussed in §|4], under the assumption that the neutron 
stars are all 1AM & and the companion masses are the central value predicted by the Pb-rri2 
relation. Of these systems, only B1855+09 has a well measured signal QRyba and Taylor 



|1991| , |Kaspi, Taylor, and Kyba 1994| ); it and the less studied J2019+2425 have the largest 
predicted signals. (It is also interesting to note that of the pulsars in Table ^, the only two 
that have been clearly shown to have "classical" interpulse emission separated from the main 
radio emission peak by ~ 180° are B1855+09 and J1804— 2718. These are two of the four 
pulsars that the -Pf>-m 2 relation suggests are observed on lines of sight that are most nearly 
equatorial, as would be expected for systems in which both magnetic poles are seen.) 

An assumed pulsar mass distribution and the -P&-m 2 relation allow us to test our assump- 
tion that the orbital inclinations are randomly distributed: i.e., cost is uniformly distributed. 
In Figure ^ we show the cumulative distribution of cos % for the thirteen pulsar-white-dwarf 
binaries discussed above, using sini values from Table ||] except for B1855+09, for which a 
better estimate is available from timing. We find that the measured values are consistent 
with uniform at the 81% level. Because the pulsar spin and orbital angular momenta are 
most likely aligned, we further conclude that there is no evidence for either alignment or 
counteralignment of pulsar magnetic fields in millisecond pulsars, contrary to some predic- 
tions (e.g., |Ruderman 199 1|) . 
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Pulsar mass (M ) 



Fig. 1. — Cumulative probabilities / dm\p{ra\' 1 f, P&) for the 13 pulsars in Table @, as de- 
scribed in the text. A vertical line is shown at 1.35M . 
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Fig. 2. — Maximum likelihood estimate of the mean and standard deviation, rh and a, of a 
gaussian neutron star mass distribution. The maximum likelihood solution is marked with 
a cross, and contours indicate 68% and 95% confidence regions. 
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1 1.2 1.4 1.6 

Minimum neutron star mass (M ) 

Fig. 3. — Maximum likelihood estimate of the minimum and maximum neutron star mass, 
to; and m u , assuming masses are uniformly distributed between the upper and lower bounds. 
The maximum likelihood solution is marked with a cross, and contours indicate 68% and 
95% confidence regions. 
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Fig. 4. — Cumulative distribution of p(f < f) (solid line, §|) and of cosi (dotted line, §|) 
for the thirteen millisecond pulsars with Pf, > 3 days and P < 10 ms, assuming a pulsar 
mass distribution ni\ = 1.35 ± O.O4M and companion mass distribution as predicted by 
the P}y-m2 relation (see text). Each distribution is consistent with uniform (dashed line, KS 
probability of 99% and 81%, respectively), as expected if orbital inclinations are randomly 
distributed and the Pfe-m 2 relation correctly predicts the companion mass distribution. 
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Fig. 5. — Neutron star masses from observations of radio pulsar systems. All error bars 
indicate central 68% confidence limits, except upper limits are one-sided 95% confidence 
limits. Five double neutron star systems are shown at the top of the diagram. In two cases, 
the average neutron star mass in a system is known with much better accuracy than the 
individual masses; these average masses are indicated with open circles. Eight neutron-star- 
white-dwarf binaries are shown in the center of the diagram, and one neutron-star-main- 
sequence-star binary is shown at bottom. Vertical lines are drawn at m — 1.35 ± O.O4M . 
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Table 1. Binary radio pulsar systems 3. 



PSR name P spin (s) P b (d) e Notes 



Double neutron star binaries 



J1518+4904 


0.040935 


8.634 


0.24948 




B1534+12 


0.037904 


0.421 


0.27368 




B1913+16 


0.05903 


0.323 


0.61713 




B2127+11C 


0.030529 


0.335 


0.68141 


1 


B2303+46 


1.066371 


12.34 


0.65837 




Neutron star /white dwarf binaries 


B0021-72E 


0.003536 


2.257 


0.000 


1 


B0021-72I 


0.003485 


0.226 


0.00 


1 


B0021-72J 


0.002101 


0.121 


0.00 


1 


J0034-0534 


0.001877 


1.589 


0.0000 




J0218+4232 


0.002323 


2.029 


0.00000 




J0437-4715 


0.005757 


5.741 


0.00000 




J0613-0200 


0.003062 


1.199 


0.00000 




B0655+64 


0.195671 


1.029 


0.00000 




J0751+1807 


0.003479 


0.263 


0.0000 




B0820+02 


0.864873 


1232.47 


0.01187 




J1012+5307 


0.005256 


0.605 


0.00000 




J1022+10 


0.016453 


7.805 


0.0001 




J1045-4509 


0.007474 


4.084 


0.00000 




B1310+18 


0.033163 


255.8 


0.002 


1 


J1455-3330 


0.007987 


76.175 


0.00017 




J1603-7202 


0.014842 


6.309 


0.0000 


2 


B1620-26 


0.011076 


191.443 


0.02531 


1,4 


J1640+2224 


0.003163 


175.461 


0.0008 




B1639+36B 


0.003528 


1.259 


0.005 


1 


J1643-1224 


0.004622 


147.017 


0.00051 




J1713+0747 


0.004570 


67.825 


0.00007 




B1718-19 


1.004037 


0.258 


0.000 


1 


B 1744- 24 A 


0.011563 


0.076 


0.0000 


1 


B1800-27 


0.334415 


406.781 


0.00051 




B1802-07 


0.023101 


2.617 


0.212 


1 
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Table 1 — Continued 



PSR name 


-fspin (s) 


Pb(d) 


e 


Notes 


J1804-2717 


0.009343 


11.129 


0.00004 


2 


B1820-11 


0.279828 


357.762 


0.79462 


3 


B1831-00 


0.520954 


1.811 


0.000 




B1855+09 


005362 


12.327 


00002 




J1910+0004 


0.003619 


0.141 


0.00 


1 


J1911-1114 


0.003626 


2.717 


0.0000 


2 


B1953+29 


0.006133 


117.349 


00033 




B1957+20 


0.001607 


0.382 


0.00000 




J2019+2425 


0.003935 


76.512 


0.00011 




J2033+17 


0.005949 


56.2 


0.00 




J2129-5721 


0.003726 


6.625 


0.0000 


2 


J2145-0750 


0.016052 


6.839 


0.00002 




J2229+2643 


0.002978 


93.016 


0.00026 




J2317+1439 


0.003445 


2.459 


0.00000 




Neutron star /main sequence 


binaries 




J0045-7319 


0.926276 


51.169 


0.808 




B1259-63 


0.047762 


1236.724 


0.86993 




Pulsar with planetary mass companions 


B1257+12 


0.006219 


66.536 


0.0182 






2nd planet 


98.223 


0.0264 





a Unless otherwise indicated, all data are from the 1995 
revision of the Princeton Pulsar Catalog (Taylor et al. 
1995). 

Note. — (1) Believed to be a globular cluster member; 
(2) Lorimer et al. 1996; (3) may be a double neutron star 



system; (4) probably a triple system (Thorsett, Arzouma- 
nian, and Taylor 1993j ). 
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Table 2. Mass estimates from the Pfc-m2 relation. 













Pulsar mass b 






Predicted: c 




A 


/ 


m 2 a 


95% 


68% 




68% 


95% 




At s 


Pulsar 


(days) 


(1O- 3 M ) 


(M & ) 


lower 


lower 


median 


upper 


upper 


sini 


M 


J0437-4715 


5.741 


1.243 


0.164 


0.106 


0.472 


1.044 


1.573 


1.971 


0.87 


3.5 


J1045-4509 


4.083 


1.765 


0.132 


0.047 


0.234 


0.557 


0.965 


1.273 


0.97 


6.2 


J1455-3330 


76.174 


6.272 


0.305 


0.098 


0.491 


1.144 


1.681 


2.058 


0.85 


5.9 


J1640+2224 


175.460 


5.907 


0.351 


0.135 


0.651 


1.495 


2.190 


2.681 


0.73 


4.6 


J1643-1224 


147.017 


0.783 


0.341 


0.553 


2.106 


4.439 


6.335 


7.674 


0.38 


1.6 


J1713+0747 


67.825 


7.896 


0.299 


0.078 


0.405 


0.960 


1.419 


1.741 


0.91 


7.6 


J1804-2718 


11.128 


4.137 


0.212 


0.073 


0.369 


0.856 


1.265 


1.558 


0.94 


6.7 


B1855+09 


12.327 


5.557 


0.219 


0.060 


0.315 


0.745 


1.101 


1.351 


0.97 


9.0 


B1953+29 


117.349 


2.417 


0.328 


0.236 


1.024 


2.251 


3.251 


3.955 


0.58 


2.8 


J2019+2425 


76.511 


10.686 


0.305 


0.062 


0.338 


0.818 


1.218 


1.500 


0.95 


10.0 


J2033+17 


56.2 


2.75 


0.290 


0.172 


0.772 


1.720 


2.489 


3.029 


0.67 


3.2 


J2129-5721 


6.625 


1.049 


0.176 


0.138 


0.591 


1.293 


1.900 


2.348 


0.80 


2.7 


J2229+2643 


93.015 


0.839 


0.315 


0.463 


1.789 


3.787 


5.409 


6.552 


0.42 


1.7 



a Central value from P^-vni relation, eqn. ^ and eqns. [l^ and 
b Median and central 68% and 95% confidence bounds. 

c Mean value of sini and Shapiro delay amplitude Ats — 2m2Tolog(l — sini), assuming a gaussian underlying 
neutron star mass distribution mi = 1.35 ± 0.04 and uniform ni2 in the range allowed by the Pb-m-i relation (§^|). 
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Table 3. Radio pulsar mass summary 



Median 68% 95% 

Star mass (Mq) central imits central limits Notcs a 



Double neutron star binaries 



J1518+4904 



pulsar 


1.56 


+0.13/ - 0.44 


+0.20/ - 1.20 


GR, RO 


companion 


1.05 


+0.45/ - 0.11 


+1.21/ - 0.14 


GR, RO 


average 


1.31 


±0.035 


±0.07 


GR 


B1534+12 










pulsar 


1.339 


±0.003 


±0.006 


GR 


companion 


1.339 


±0.003 


±0.006 


GR 


B1913+16 










pulsar 


1.4411 


±0.00035 


±0.0007 


GR 


companion 


1.3874 


±0.00035 


±0.0007 


GR 


B2127+11C 










pulsar 


1.349 


±0.040 


±0.080 


GR 


companion 


1.363 


±0.040 


±0.080 


GR 


average 


1.3561 


±0.0003 


±0.0006 


GR 


B2303+46 










pulsar 


1.30 


+0.13/ - 0.46 


+0.18/ - 1.08 


GR, RO 


companion 


1.34 


+0.47/ - 0.13 


+1.08/ - 0.15 


GR, RO 


average 


1.32 


±0.025 


±0.05 


GR 


Neutron star/white dwarf binaries 


J0437-4715 






< 1.51 


x, P b m 2 


J1012+5307 


1.7 


±0.5 


±1.0 


Opt 


J1045-4509 






< 1.48 


P b m 2 


J1713+0747 


1.45 


±0.31 


±0.62 


P b m 2 , GR 




1.34 


±0.20 


±0.40 


P b m 2 , GR, Opt 


B1802-07 




< 1.39 


< 1.45 


GR 




1.26 


+0.08/ - 0.17 


+0.15/ - 0.67 


GR, P b m 2 


J1804-2718 






< 1.73 


P b m 2 


B1855+09 


1.41 


±0.10 


±0.20 


GR 


J2019+2425 






< 1.68 


P b m 2 


Neutron star/main sequence binaries 


J0045-7319 


1.58 


±0.34 


±0.68 


Opt, MS 



a Assumptions made in mass estimate: (GR) general relativistic binary model; 
(RO) random orbital orientation, or inclination angle uniform in cosi; (P b m 2 ) core 
mass orbital period relation; (x) proper-motion-induced change in the projected semi- 
major axis; (Opt) optical companion observations; (MS) main sequence stellar model. 



